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We preform a multi-orbital analysis on the novel superconductivity in Na x Co02 ■ J/H2O. 
We construct a three-orbital model which reproduces the band structure expected from the 
LDA calculation. The effective interaction leading to the pairing is estimated by means of the 
perturbation theory. It is shown that the spin triplet superconductivity is stabilized in the wide 
parameter region. This is basically owing to the ferromagnetic character of spin fluctuation. The 
p-wave and /-wave superconductivity are nearly degenerate. The former is realized when the 
Hund's rule coupling is large, and vice versa. In a small part of the parameter space, the d-wave 
superconductivity is also stabilized. We point out that the orbital degeneracy plays an essential 
role for these results. In particular, the momentum dependence of wave function of quasi- 
particles is quite important. The nearly degeneracy of p-wave and /-wave superconductivity is 
explained by analysing the microscopic character of each Fermi surface. We discuss the validity 
of some reduced models. While the single-orbital Hubbard model reproducing the Fermi surface 
is qualitatively inappropriate, we find an effective two-orbital model appropriate for studying 
the superconductivity. We investigate the vertex corrections higher than the third order on the 
basis of the two-orbital model. It is shown that the vertex correction induces the screening 
effect but does not affect qualitative results. 

KEYWORDS: Na x Co02 -yHbO; unconventional superconductivity; ferromagnetic spin fluctuation; multi- 
orbital analysis 



1. Introduction 

Since the discovery of High-T c superconductivity 1 -' 
and heavy fermion superconductors, 2 ' the mechanism 
of superconductivity induced by electron correlation has 
been one of the central issues in the condensed matter 
physics. In this study, recently discovered superconduc- 
tor Na x Co02 • J/H2O is analyzed in details. 

Immediately after the discovery of superconductivity 
in water-intercalated Cobalt oxides Na x Co02 • 2/H2O, 3 ) 
both experimental 4-20 ) and theoretical 21-32 ) studies 
have been performed extensively. While some contro- 
versial results exist, many experimental evidences for 
the non-s-wave superconductivity 20 -' has been reported 
by NMR 10-13 ' and specific heat measurements. 17-19 - 1 
The characteristic behaviors in strongly correlated elec- 
tron systems have been observed in the non-water- 
intercalated compounds. 7-9 ' 33 ) The existence of the mag- 
netic phase 4-6 ) in Na x Co02 with x ~ 0.75 also indicates 
an importance of electron correlation. These compounds 
have a layered structure like cuprate 1 ) and ruthenate, 34 ) 
and the two-dimensionality is enhanced by the water- 
intercalation. These circumstantial evidences indicate 
that Na x Co02 -yH^O is an unconventional superconduc- 
tor induced by the electron correlation. 

The theoretical interests are turned on also by the sym- 
metry of crystal structure. In contrast to the square lat- 
tice in cuprates and ruthenates, the layer is constructed 
from the triangular lattice of Co ions. Then, a novel sym- 
metry of Cooper pairing is possible in principle. The d- 
wave superconductivity in cuprate superconductors and 
p-wave superconductivity in ruthenates have been estab- 
lished before. In addition to them, the spin triplet /-wave 
superconductivity and spin singlet i-wave one are possi- 
ble from the analysis of pairing symmetry (see Table. I). 



The effect of frustration, which is characteristic in the 
spin system on the triangular lattice, has also attracted 
much attention. The RVB theory has been applied to the 
triangular lattice 23-26 ) and basically concluded the spin 
singlet d-wave superconductivity. Then, d x 2_ y 2± ic£ xy - 
wave symmetry is expected below T c owing to the six- 
fold symmetry of triangular lattice. However, the time- 
reversal symmetry breaking has not been observed until 
now. 14 ' Some authors have pointed out the frustration 
of charge ordering for the electron filling n = 4/3, 25 ) and 
the /-wave superconductivity due to the charge fluctua- 
tion has been discussed. 29 ' 32 ) 

Another interesting property of Na x Co02 -yl^O is the 
orbital degeneracy. The conduction band of this material 
mainly consists of three t2g-orbitals in Co ions which hy- 
bridize with 02p-orbitals. Thus far, most of theoretical 
studies on the superconductivity have been performed 
on the basis of the single-orbital model. These investiga- 
tions have successfully achieved microscopic understand- 
ings on the cuprate, organic and ruthenate superconduc- 
tors. 35 ) However, we consider that the theoretical anal- 
ysis including the orbital degeneracy is highly desired in 
order to understand a variety of superconductors includ- 
ing Na x Co02 • J/H2O and heavy fermion compounds. The 
superconductivity in ri-electron systems provides a favor- 
able subject for the theoretical development along this 
line, because a simple electronic structure is expected 
compared to heavy fermion superconductors. Although 
Sr2Ru04 has been a precious compound in this sense, 
then the orbital degree of freedom is not important for 
the basic mechanism of superconductivity 36,37 ) In this 
study, we show that the orbital degeneracy plays an es- 
sential role in Na x Co02 -j/HaO in contrast to the ruthen- 
ate superconductor. We conclude that Na x Co02 ■ J/H2O 
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is a typical multi-orbital superconductor in this sense. 

We adopt a perturbative method for the unconven- 
tional superconductivity, 35 ^ which is a systematic ap- 
proach for the electron correlation. Note that the spin 
fluctuation theory 38 ) which is widely used for supercon- 
ductivity is microscopically formulated in this method. 
It is expected that this approach is reliable from weak 
to intermediate coupling region. Before the discovery 
of Na x Co0 2 • 2/H 2 0, this method has been applied to 
the single-orbital triangular lattice model. Then, the d- 
wave, 39 ) /-wave 40 ) and p-wave superconductivity 41 ) have 
been obtained. Some authors have applied this calcula- 
tion to Na x Co0 2 • ?yH 2 0, and reported the spin singlet 
d- or i-wave superconductivity, 27 ) spin triplet /-wave su- 
perconductivity 30 ) and nearly degeneracy between d- and 
/-wave superconductivity. 31 ) We consider that this puz- 
zling problem should be resolved by the multi-orbital 
analysis involving the microscopic aspects of electronic 
structure. 

In this paper, we analyze a multi-orbital Hubbard 
model constructed from three Co i 2g -orbitals. This model 
appropriately reproduces the electronic structure ob- 
tained in the LDA calculation. 42 ' 43 ) The wave function 
of quasi-particles, which is neglected in the single-orbital 
Hubbard models, is appropriately taken into account 
in this multi-orbital model. We show that the momen- 
tum dependence of this wave function plays an essen- 
tial role for the mechanism of superconductivity. We de- 
termine the most stable superconducting state with use 
of the perturbation theory. According to the results of 
second order perturbation (SOP), third order perturba- 
tion (TOP) and renormalized third order perturbation 
(RTOP) theories, it is concluded that the spin triplet p- 
wave or /-wave superconductivity is stable in the wide 
region of parameter space. The pairing interaction is 
closely related to the ferromagnetic character of spin sus- 
ceptibility, although the pairing interaction is not simply 
described by the spin susceptibility like in the single- 
orbital model. 35 ) While the momentum dependence of 
spin susceptibility is usually not remarkable in the frus- 
trating system, the ferromagnetic character clearly ap- 
pears in the present case owing to the orbital degree of 
freedom. 

From a comparison with single-orbital Hubbard mod- 
els, the important roles of orbital degeneracy are illu- 
minated in §4.1. Alternatively, we propose a reduced 
two-orbital model including the e g -doublet in §4.2. It is 
shown that results for the superconductivity is appropri- 
ately reproduced in this simplified model. On the basis of 
the two-orbital model, we investigate the roles of vertex 
correction terms in §5. Then, we show that the vertex 
correction term, which significantly enhances the spin 
triplet pairing in Sr 2 Ru04, 36 ) is not important in case 
of Na x Co02 • J/H 2 0. Thus, the superconducting insta- 
bility is basically described within the SOP. Therefore, 
we first explain in details the results of SOP in §3, and 
discuss the reduced models in §4 and the role of vertex 
corrections in §5. 



2. Multi-orbital model 

First, we construct a multi-orbital model for 
Na x Co0 2 • yH 2 44 ) We consider a two-dimensional 
model which represents the Co ions on the trian- 
gular lattice. Note that the superconductivity occurs 
when the two-dimensionality is enhanced by the water- 
intercalation. We also note that the conduction band 
mainly consists of Co i 2g -orbitals. 42,43 ) Co ion is enclosed 
by an octahedron of oxygens and nearest neighbor Co 
ions share the edge of the octahedron. We describe the 
dispersion relation by using a tight-binding model and 
adopt a multi-orbital Hubbard Hamiltonian written as, 

H 3 = H Q + H It (1) 

Ha= ^2^2 ta , b ^<3 C la,s C 3,b,s, (2) 
i,j,s a.b 

H l = U ^2 ^2 ni,a,tn it a,l + U' ^ ^2 n i^ n i,b 
i a i a>b 

~Jh "^2 "^2(2Si,aSi,b + X"i,a«-a) 



i a>b 



+J J2J2 C la,l C la,f c iM c iM- 



(3) 



The first term H a is a tight-binding Hamiltonian where 
ta,b,i,j are hopping matrix elements. Here, the indices i 
and j denote the sites in the real space and indices a 
and b denote the orbitals. We assign the d xy -, d yz - and 
d xz -orbitals to a = 1, a = 2 and a = 3, respectively. 
The largest matrix element is the inter-orbital hopping 
through 02p-orbitals, which are ti,2,»,j for j = i±(a+b), 
h,3,i,j for j = «±a and t\^,i,j f° r 3 = i±b. We choose the 
lattice constant as a unit length and denote the unit vec- 
tors as a=(\/3/2, —1/2) and b=(0, 1) which are the basis 
of the triangular lattice. If we assume only the largest 
matrix elements, the system is regarded to be a superpo- 
sition of the kagome lattice. 22 ) However, the long range 
hopping through the 02p-orbitals and direct hopping be- 
tween Co ions are necessary to reproduce the Fermi sur- 
face obtained in the LDA calculation. 

We take account of the matrix elements within third- 
nearest-neighbor sites according to the symmetry of or- 
bitals and lattice. They are described by nine parame- 
ters from t\ to tg. The non- interacting Hamiltonian is 
described in the matrix representation, 



(4) 



en(fe) £i2 (fc) £13 (*0 
H(k) = ( £ 2 i(fe) s 22 (k) e 23 (fc) I , (5) 

£3l(fe) £32(fc) £ 3 3(fc) 

where ct = (ct , ct , ct ) is a vector representa- 

k,s v fe,l,s' fe,2,s' k,3,s' * 

tion of the Fourier transformed creation operators with 
spin s. The matrix elements are obtained as, 

£n(fe) = 2£i cosfci + 2i 2 (cosfc 2 + cosk 3 ) 

+2i4(cos(fci — k 3 ) + cos(fci — fc 2 )) + 2t 5 cos2fci, (6) 

£22 (k) — 2ti cosfc 2 + 2£ 2 (cosfci + cosfcs) 
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-r-2i 4 (cos(fci - k 2 ) + cos(fc 2 - fc 3 )) + 2t 5 cos2fc 2 , (7) 

£33 (fe) = 2ti cosfc 3 + 2i 2 (cosfci + cosfc 2 ) 

+2t 4 (cos(fci - fc 3 ) + cos(fc 2 - fc 3 )) + 2t 5 cos2fc 3 , (8) 

£i2(fe) = 2i 3 cosfc 3 + 2i 6 cos2fc 3 + 2tr cos(£:i — fc 3 ) 

+2t$ cos(fc 2 - fc 3 ) + tg cos(fci - fc 2 ) - e c /3, (9) 

£ i3(fc) = 2t 3 cosfc 2 + 2i 6 cos2fc 2 + 2t 7 cos(k 2 — fc 3 ) 

+2t 8 cos(fc! — fc 2 ) + tg cos(fci — fc 3 ) — e c /3, (10) 

£2 3 (fc) = 2i 3 cosfci + 2t 6 cos2fci + 2t?cos(ki — fc 2 ) 

+2i 8 cos(A;i - fc 3 ) + t 9 cos(fc 2 - h) - e c /3, (11) 

where fci = V3/2k K — l/2fc y , fc 2 = fc y and fc 3 = — fci — fc 2 . 
The parameter e c represents the crystal field splitting of 
i 2g -orbitals arising from the distortion of octahedron. A 
typical dispersion relation and Fermi surface are shown 
in Fig. 1. There is a hole pocket enclosing the L-point 
and six hole pockets near the K-points, which arc con- 
sistent with LDA calculations. 42 ' 43 ^ We choose the unit 
of energy as t 3 = 1 throughout this paper. 

Although e c seems to be small, it is useful to use 
a non-degenerate ai g -orbital and doubly-degenerate e g - 
orbitals. They are defined from the three i 2g -orbitals as 

>= -iflxz > -|yz >), 



(a) 



|e g ,l 



V2' 



,2>= 



|ffli g >= 



7! (2|xy> 



xz > 



|yz >), 



-^=(|xy > +|xz > +|yz >). 



(12) 



(13) 



(14) 



The wave function of ai g -orbital spreads along the c-axis, 
and those of e g -orbitals spread along the two-dimensional 
plane. We will show later that this representation is ap- 
propriate for understanding the mechanism of supercon- 
ductivity (§4.2). 

The hole pocket around the r-point in Fig. 1(a) mainly 
consists of the a lg -orbital and the six hole pockets near 
the K-points mainly consist of the e g -orbitals. Thus, 
we denote these Fermi surfaces as a lg -Fermi surface 
and e g -Fermi surface, respectively. This nature of the 
Fermi surface is consistent with LDA calculations. 42,43 ' 
Note that recent ARPES measurements 45,46 ' for non- 
superconducting Na^Co0 2 observed the ai g -Fermi sur- 
face, but the e g -Fermi surface has not been found. Fermi 
surface of water-intercalated Na x Co0 2 is not clear at 
present. Moreover, the valence of Co ion in superconduct- 
ing materials is also under debate. 47 ' Therefore, we in- 
vestigate a wide region in the parameter space and study 
the possible pairing instability. It is one of the goals of 
this paper to study the relation between the electronic 
state and superconductivity. It will be shown that the 
superconductivity is hard to be stabilized when e g -Fermi 
surface vanishes. 

The second term H\ describes the short range Coulomb 
interactions which include the intra-orbital repulsion U , 
inter-orbital repulsion U', Hund's rule coupling Jh and 
pair hopping term J. The relations U = U' + Jh + J and 
Jh = J are satisfied in a simple estimation. Under these 
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Fig. 1. (a) Fermi surfaces and (b) dispersion relation ob- 
tained from the tight-binding Hamiltonian. The dashed 
line in (a) shows the first Brillouin zone. The pa- 
rameters are chosen to be (ti , ti, *3, t4, ts, te, <7, tg, tg) = 
(0.08, 0.16, 1, 0.24, -0.16, -0.04, 0.16, 0.16, -0.2). 



conditions, the interaction term Hi is invariant for the 
local unitary transformation between orbitals which will 
be used later. If these relations are violated, the symme- 
try of triangular lattice is artificially broken. Therefore, 
we impose these relations through this paper. Although 
possible roles of the long range Coulomb interaction have 
been investigated, 23 ' 29 ' 32 ' we concentrate on the short 
range interaction in this paper. 

Note that previous studies based on a perturbative 
method for cuprates, organics and ruthenate have suc- 
ceeded in identifying the dominant scattering process 
leading to the superconductivity. 35 ' This theory is com- 
plementary to the fluctuation theory which is represented 
by a random phase approximation (RPA) or fluctuation 
exchange approximation (FLEX). Generally speaking, 
the fluctuation theory will be appropriate in the vicinity 
of the magnetic or other instabilities, because the critical 
enhancement of the fluctuation is taken into account. On 
the other hand, the perturbation theory is more appro- 
priate when the critical enhancement of any particular 
fluctuation is absent, because all terms in the same order 
are taken into account without any prejudice. We per- 
form the second order perturbation as well as the third 
order perturbation in this paper. The results of FLEX 
study will be published elsewhere, 44 ' where qualitatively 
consistent results are obtained. 
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3. Second Order Perturbation 

3. 1 Details of calculation and classification of pairing 
symmetry 

In this section, we investigate the superconducting in- 
stability by using the Eliashberg equation within the 
second order perturbation (SOP). The basic procedure 
has been explained in literatures 35 ^ and the extension to 
multi-orbital model is straightforward. The Eliashberg 
equation is described by the Green function and the ef- 
fective interaction. The latter is represented by an irre- 
ducible four point vertex in the particle-particle channel 
(Fig. 2(a)). The second order terms in the effective inter- 
action are diagrammatically represented by Figs. 2(b-e). 
In case of the single-orbital Hubbard model, this term 
is simply expressed as V(k, k') = U 2 xo(k — k') for spin 
singlet pairing and V(k,k') = —U 2 xo(k — k') for spin 
triplet pairing, respectively, with a bare spin susceptibil- 
ity Xo{k — k'). However, in the multi-orbital model, the 
four point vertex has indices of orbitals as V a b c d(k,k') 
(see Fig. 2(a)), which is calculated from the possible com- 
bination of Coulomb interactions and Green functions. 



(a) 
a 



(b) 



(c) 



b->- 





-) — ; — >- 



(d) 



(e) 



) * ■ ■ ) 



-> : >- 



Fig. 2. (a) Diagrammatic representation of the effective interac- 
tion leading to the superconductivity, (b-e) The second order 
terms with respect to the Coulomb interactions (dashed lines). 
The solid line denotes the Green function having the indices of 
spin and orbital. 



In order to make the following discussions clear, we in- 
troduce a unitary transformation U(k) — (uij(k)) which 
diagonalizes H(k), namely 



(j\k)H{k)U{k) 




(15) 



Here, we choose -Ei(fc) < E 2 (k) < E 3 {k). With use of 
these matrix elements, the matrix form of Green func- 
tion characterized by orbitals G(k) = (iw„l — -ff(fc)) -1 
is described as, 



Gij(k) = ^u ta (k)u ja (k)G a {k), 



(16) 



1 



where G a (k) = 

In the following, we denote the energy band described 
by the dispersion relation E 3 {k) as 7-band. As we have 
shown in Fig. 1, the 7-band crosses the Fermi level, and 
the others are below the Fermi level. Therefore, the su- 



perconducting transition is induced by the Cooper pair- 
ing in the 7-band. In this case, the Eliashberg equation 
is written in terms of an effective interaction within the 
7-band, 

A c A(fc) = V ^ k')\G 3 (k')\ 2 A(k>), (17) 
k> 

with 

V(k,k') = Y u a3 (k)u b3 (-k)V abc d{k,k')u c3 (k')u d3 (-k'). 

(18) 



abed 



The Eliashberg equation (eq. (17)) is regarded to be 
an eigenvalue equation and A e represents the maximum 
eigenvalue. The superconducting transition temperature 
is determined by the criterion A c = 1. 

Here, we have ignored the normal self-energy which 
is important for a quantitative estimation of T c . How- 
ever, qualitative nature of the superconductivity, such 
as the pairing symmetry and the pairing mechanism, is 
not affected in many cases including cuprates, ruthen- 
ates and organics. 35 -* This is highly expected in case of 
Na x Co02 • 2/H2O, unless the electronic structure is sig- 
nificantly affected by the normal self-energy. We will 
show that the volume of e g -Fermi surface , which will 
be denoted as n e below, is an important parameter for 
the pairing symmetry. Therefore, it is possible that the 
pairing symmetry is affected by the normal self-energy 
through the modification of n c . It is, however, expected 
that the following results are still valid even in this case 
by regarding the n c modified by the normal self-energy 
as a relevant parameter. 



Ai 


s-wave 


1 


E2 


d-wave 


sin -^fc x sin ^k y 


A 2 


i-wavc 


sin ^5 fe x sin ifc y + sin ^fc x sin |fc y — sin \/3fcx sm 2fcy 


Ei 


p-wave 


sin fcx cos ^ k Y 


Bi 


/i-wave 


1 JT 1 — " 

sin i fcy (cos fc x — cos i fc y ) 


B 2 


/2-wave 


sin fc x (cos fc x - cos |fc y ) 



Table I. Classification of the pairing symmetry in the triangular 
lattice. The first column shows the irreducible representations 
of Dg group. The second column shows the notation adopted 
in this paper. The s-wave, p-wave, etc are the counterparts of 
the isotropic system. The third column shows the typical wave 
function of Cooper pairs. 



Before showing the results, it is necessary to classify 
the pairing symmetry. The symmetry of Cooper pairs is 
classified into s-, p- 7 d-wave etc. in case of an isotropic 
system like 3 He. For metals, the Cooper pairing is classi- 
fied into the finite species according to the symmetry of 
crystals. 48 ^ We show the classification in case of the tri- 
angular lattice in Table I. We denote "s-wave" , "d-wave" 
etc. in analogy with the isotropic case. While the s-, al- 
and i-wave are spin singlet pairings, the p-, f\- and / 2 - 
wave are spin triplet pairings. Note that there remains 
two-fold degeneracy in the p- and cZ-wave states, namely 
p x - andpy-wave, d xy - and d x 2_ y 2-wave, respectively. The 
time-reversal-symmetry-breaking is expected below T c in 
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the d-wave state, as discussed in the RVB theory. 23 ~ 26 ) 
On the contrary, time-reversal-symmetry is not necessar- 
ily broken in the p-wave case because there is an internal 
degree of freedom representing the direction of S = 1, as 
discussed in Sr 2 Ru04. 37 ^ 

The eigenvalues of the Eliashberg equation, eq. (17) 
are classified according to the symmetry of Cooper 
pairs. The pairing symmetry corresponding to the largest 
eigenvalue is stabilized below T c . Hereafter, we ignore the 
possibility of s-wave pairing because the strong on-site 
repulsion will destabilize even the extended s-wave pair- 
ing. When the symmetry of crystal is lowered, some can- 
didates in Table I are classified into the same irreducible 
representation. For example, the d xy -wave and s-wave 
symmetries are included in the same representation for 
the anisotropic triangular lattice. 49 ' 50 ) However, we can 
ignore this possibility in the isotropic triangular lattice. 

3.2 Phase diagram of three- orbital model 

In order to search possible pairing symmetries in a 
phase diagram, we introduce two controlling parameters, 
a and n c . Among the hopping matrix elements in eqs. (6- 
11), the largest one, namely £3 is fixed to 1 but the other 
matrix elements are chosen to be 

(*1>*2>*4)*5>*6>*7>*8>*9) = 

o(0.1, 0.2, 0.3, -0.2, -0.05, 0.2, 0.2, -0.25). (19) 

We choose this parameter set so that the dispersion re- 
lation obtained in the LDA calculation 42 ' 43 ) is appro- 
priately reproduced when a <~ 1. In case of a = 0, the 
system is regarded to be a superposition of kagome lat- 
tice, 22 ) but we have to choose a > 0.6 in order to obtain 
a realistic Fermi surface. Thus, the parameter a indicates 
a deviation from the kagome lattice. Although there are 
many choices of controlling the minor matrix elements, 
we have confirmed that the following results are qualita- 
tively independent of the choice. 

As another controlling parameter, we use the hole 
number n c in the e g -Fermi surface, which can be altered 
by adjusting the crystal field splitting e c . When we de- 
crease e c , the energy of e g -orbitals is lowered and thus 
n e decreases. We have confirmed that the value n c is es- 
sential rather than the total electron number n for the 
following results which are almost independent of the 
way to alter n e . Note that the total electron number is 
fixed as n = 5.33 throughout this paper. 

We divide the first Brillouin zone into 128 x 128 lat- 
tice points and take 512 Matsubara frequencies. We have 
confirmed that the following results do not depend on 
the numerical details, qualitatively. In the following, the 
temperature is fixed to be T = 0.01 unless we mention 
explicitly. It will be shown in Fig. 5 that the stable pair- 
ing symmetry is almost independent of the temperature. 
We fix U = 5 and change the value of Jh = J. Under the 
reasonable conditions U = U' + 2 Jh and U' — Jh > 0, 
Jh = U/3 is the maximum value of the Hund's rule cou- 
pling. 

Figure. 3 shows the most stable pairing symmetry in 
the phase diagram of a and n e for two values of the inter- 
action strength. As shown in Fig. 3(a), the spin triplet 




0.1 0.2 0.3 0.4 0.5 0.6 



n 

e 




0.1 0.2 0.3 0.4 0.5 0.6 

n 

e 



Fig. 3. Phase diagram for (a)CI' = J H = J = U/3 and (b)U' = 
U/2 and Jjj = J = U/4. The horizontal and vertical axes are de- 
scribed in the text. The solid line is the phase boundary obtained 
by the interpolation. 



p-wave superconductivity is stabilized in the wide region 
of parameter space when Jh = U /3. The /i-wave super- 
conductivity is also stabilized when e g -Fermi surface is 
very small or very large. For the values of n c expected in 
the LDA calculation, namely n c = 0.1 ~ 0.3, we obtain 
the p-wave superconductivity independent of the value of 
a. When the value of Hund's rule coupling is decreased 
(Fig. 3(b)), the /i-wave superconductivity becomes more 
stable. We see that in both cases the spin triplet super- 
conductivity is stable. 

By definition, the e g -Fermi surface vanishes in case 
of n c = 0. Then, it is difficult to determine the pair- 
ing state since the tendency to superconductivity is very 
weak independent of the pairing symmetry. On the other 
hand, the superconductivity is not significantly affected 
by the disappearance of ai g -Fermi surface which occurs 
at n c = 0.67. 

In order to make the situation clearer, we show the 
eigenvalues of Fliashberg equation for each pairing sym- 
metry in Fig. 4. It is shown that the p- and /i-wave 
superconductivity have nearly degenerate eigenvalues in 
a wide parameter range. If we assume the weak crys- 
tal field splitting e c ~ 0, we obtain n c ~ 0.3 which is 
consistent with LDA calculation. The eigenvalue in the 
/i-wave symmetry shows a minimum around this value. 
As a result, the p-wave superconductivity is stable in this 
region. As the Hund's rule coupling decreases, eigenval- 
ues of both p- and /i-wave symmetries increase, but that 
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n 

e 

Fig. 4. n c -depcndence of eigenvalues of Eliashberg equation. We 
choose T = 0.01, a = 0.8 and (a) 17' = Jh = J = U/3 or (b)C7' = 
[7/2 and J H = J = 17/4. 



of the /i-wave symmetry increases more rapidly (See also 
Fig. 7). Note that the eigenvalues for the (i-wave, i-wave 
and /2-wave states are very small compared to the p- and 
/i-wave states. As is shown later, the d-wave state is sta- 
bilized when Hund's rule coupling is very small (Figs. 7 
and 8). 
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Fig. 5. Temperature dependence of eigenvalues of Eliashberg 
equation. We choose a = 0.8, n c = 0.238, U' = U/2 and 
J H = J = [7/4. 



We see that A c is still less than 1 at T = 0.01 (Fig. 4). 



Therefore, the pairing instability occurs at lower temper- 
ature. Fig. 5 shows the temperature dependence of A c at 
Jh = U/4, a = 0.8 and n c = 0.238 where the maximum 
eigenvalue is A c <~ 0.7 at T = 0.01. Then, we obtain 
A c = 1 at T c = 0.0037 for the p-wave symmetry. If we 
assume t 3 = 200meV so that the total band width is 
W = 1.8eV, T c = 0.0037 corresponds to T c = 8K consis- 
tent with experimental value. Furthermore, Fig. 5 clearly 
shows that most stable pairing symmetry is almost in- 
dependent of temperature. This means that the phase 
diagram obtained at T = 0.01 is very accurate. 

Another interesting result in Fig. 4 is that the max- 
imum eigenvalue does not significantly depend on n c . 
Even if the size of e g -Fermi surface is remarkably re- 
duced, the instability of superconductivity is not sup- 
pressed unless the e g -Fermi surface vanishes. This is 
mainly because the DOS of e g -Fermi surface little de- 
pends on the value of n e . This is one of the characteristics 
of the two-dimensional system in the low density region. 
Note that the number of hole included in each hole pocket 
is very small as n c /6 ~ 0.05. Then, an analogy with the 
isotropic system like 3 He is partly justified. This picture 
is important for the pairing mechanism as we will ex- 
plain in §3.3. The n c -dependence of T c can be measured 
by varying the Na-content of Na x Co02 • J/H 2 0. However, 
experimental results seems to be controversial. 51,52 ' 

The eigenvalue rapidly decreases when e g -Fermi sur- 
face vanishes. This result indicates that the e g -Fermi 
surface plays an essential role for the superconductiv- 
ity. This implication will be clearly confirmed in §4.2. 
Although the eigenvalues are very small, the cZ-wave sym- 
metry seems to be most stable at n c = 0. Then, the topol- 
ogy of Fermi surface is equivalent to the simple triangu- 
lar lattice including only the nearest neighbor hopping. 
In this sense, our result at n c — is qualitatively consis- 
tent with the RVB theory based on the t-J model in the 
triangular lattice, which shows the d x 2_ y 2 ±ici xy -wave su- 
perconductivity. 23-26 ) However, the used parameters are 
quite different. The t-J model assumes U/t > 8, while 
U/t = 5 in this paper. In the intermediate coupling re- 
gion, the momentum dependence arising from the vertex 
correction is probably important when the SOP gives 
very small A c . 35 ) In case of the simple triangular lat- 
tice, the lowest order vertex correction favors the p-wave 
state. 41 ) It should be stressed that the SOP gives much 
larger value of A when e g -Fermi surface exists, as shown 
in Fig. 4. 

Fig. 6 shows the a-dependence of eigenvalues. It is 
shown that the eigenvalue monotonically increases with 
decreasing a. This variation is basically owing to the in- 
crease of the DOS. In case of a = 0.5, almost flat band is 
realized around the e g -Fermi surface. Therefore, a steep 
increase of the eigenvalue leading to the remarkable en- 
hancement of T c occurs toward a — 0.5. We note that 
most important parameter for the appearance of flat 
band is the next nearest neighbor hoppings. Although by 
changing the parameter a, the nearest and third nearest 
neighbor hoppings vary simultaneously, these parameters 
play only quantitative roles. From Figs. 3-6, we see that 
the variable a is important for the value of T c , while the 
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from n c = 0.33 at a = 0.6 to n e = 0.31 at a = 1. 
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line is the phase boundary obtained by the interpolation. 



variable n c plays an essential role for determining the 
pairing symmetry. 
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Fig. 7. JH-dcpcndcncc of eigenvalues of Eliashbcrg equation. The 
parameters are chosen to be a = 0.8 and n c = 0.31 



Before closing this subsection, let us discuss the pos- 
sibility of c?-wave superconductivity in case of the small 
Hund's rule coupling. Fig. 7 shows the Jn-dependence of 
eigenvalues for each pairing symmetry. It is shown that 
all eigenvalues increase with the decrease of Hund's rule 
coupling. Among them, the eigenvalue in the d-wave sym- 
metry increases most rapidly and the d-wave supercon- 
ductivity is stabilized for Jh < U/12. The phase diagram 
in the Ju-n c plane is shown in Fig. 8. We see that the d- 
wave superconductivity is more stable when n e is small. 

This stability of the ci-wave pairing is basically ow- 
ing to the large value of U' which is comparable to U . 
The inter-orbital repulsion V couples to the charge and 
orbital excitations which contribute to the effective inter- 
actions equivalently in the singlet and triplet channels. 
Therefore, the difference between singlet and triplet su- 
perconductivity is reduced when U' is large. In other 
words, the Hund's rule coupling favors the spin triplet 
superconductivity, although the value of T c is reduced. 
However, we expect that the d-wave superconductivity is 
less stable if we include the higher order terms because 



higher order terms significantly enhance the spin exci- 
tation rather than the orbital and charge excitation. In 
other words, the role of U' will be reduced in the higher 
order theory. This is confirmed by the FLEX calcula- 
tion. 44 ) 

3. 3 Basic mechanism of superconductivity 

In order to clarify the basic mechanism of supercon- 
ductivity, we study the momentum dependence of effec- 
tive interaction V(k, k') in the spin triplet channel. Fig- 
ure 9 shows the fc'-dependence of V(k, k') with k being 
fixed at the momentum shown by an arrow at which the 
order parameter in the p-wave symmetry takes maximum 
value. It is apparent that there is a strong attractive in- 
teraction between momenta included in the same hole 
pocket Fermi surface. This is the reason why the spin 
triplet superconductivity is favored. We can show that in 
case of Jh = U/3, the effective interaction in the singlet 
channel has opposite sign to that in the triplet channel. 
This strong repulsive interaction remarkably suppresses 
the spin singlet superconductivity. 

The microscopic origin of this momentum dependence 
can be understood as follows. First, we point out the 
ferromagnetic character of spin fluctuation. Fig. 10(a) 
shows the spin susceptibility which is estimated by the 
Kubo formula within the bubble diagram. It is clearly 
that the spin susceptibility has a trapezoidal peak around 
q = 0. Note that the ferromagnetic spin fluctuation has 
been expected in the LDA calculation 42 ) and observed by 
the NMR measurement. 13 ) Owing to the ferromagnetic 
character of spin susceptibility, the attractive interaction 
in the same hole pocket is very strong and favors the spin 
triplet superconductivity. 

The ferromagnetic spin fluctuation is basically comes 
from the e g -Fermi surface. Each hole pocket gives rise 
to the ferromagnetic spin fluctuation like in the two- 
dimensional electron gas, which has a susceptibility 
with the trapezoidal structure. Actually, as shown in 
Fig. 10(a), when we increase the size of hole pockets by 
changing n c , the width of the trapezoidal peak around 
the T point increases. 

Next, we illuminate the essential roles of the orbital 
degree of freedom. First, we point out that the ferro- 
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Fig. 9. Contour plot of the effective interaction V(k, k 1 ). The ini- 
tial momentum k is shown in the figure. The horizontal and ver- 
tical axis show fc^ and k' y , respectively. Matsubara frequency is 
fixed to the lowest value oj n = uj' n = irT. The Fermi surface is si- 
multaneously described by the thin solid line. The parameters are 
chosen to be n c = 0.31, a = 0.8, U' = U/2 and J H = J = U/A. 



magnetic spin fluctuation is indeed induced by the or- 
bital degree of freedom. In the multi-orbital model, the 
spin susceptibility is determined by the dispersion re- 
lation and the structure factor arising from the orbital 
degree of freedom. If we neglect the momentum depen- 
dence of structure factor as was done in the previous 
studies, 31,53 ) we obtain two peaks of spin susceptibility 
which are quit different from ours. One is located around 
the M point and the other is slightly removed from the 
T point. However, we obtain the trapezoidal peak cen- 
tered at the r point by appropriately taking account of 
the structure factor. Thus, the frustration inherent in 
the triangular lattice is removed by the orbital degree of 
freedom which gives rise to the ferromagnetic spin fluc- 
tuation. 

Second, we point out that the roles of the orbital de- 
gree of freedom can be understood by considering the 
momentum dependence of the wave function which is 
expressed by the unitary matrix U(k) in eq. (15). This 
wave function indicates the orbital character of quasi- 
particles (see also §4). The structure factor of spin dis- 
cussed above is also obtained by this wave function. 
Furthermore, the effective interaction V(k, k') has an- 
other distinct property arising from this momentum de- 
pendence. As we have mentioned before, the e g -Fermi 
surface mainly consists of the e g -doublet whose wave 
function is shown in eqs. (12) and (13). Furthermore, 
we find that the six hole pockets are divided into three 
pairs as is shown in Fig. 10(b). For example, more than 
90% of the weight of wave function in the Fermi surface 
"A" originates from the orbital |e g , 1 >, while the other 
two pairs are dominated by respective linear combina- 
tions of |e g , 1 > and |e g ,2 >. It is generally expected 
that the electron correlation between the same orbitals 
is stronger than that between the different orbitals. Ac- 
tually, the effective interaction between different pairs 
"A" , "B" and "C" is significantly smaller than those be- 




Fig. 10. (a) Momentum dependence of the static spin suscepti- 
bility at a = 0.8. (b) Schematic figure for the classification of 
hole pockets. 



tween the same pairs, as shown in Fig. 9. This is the 
reason why the p- and /i-wave superconductivities are 
stabilized with nearly degenerate eigenvalues as shown 
in Fig. 4. Which is more stable between p- and /i-wave 
states depends on the coupling between different pairs 
of hole pockets, which is generally small as explained 
above. Note that if we apply the phenomenological the- 
ory on the ferromagnetic spin-fluctuation-induced super- 
conductivity to Na x Co02 • J/H 2 0, the /i-wave supercon- 
ductivity is much more stable rather than the p-w&ve 
superconductivity. The single band model leading to the 
ferromagnetic spin fluctuation 30 ) also concludes the f\- 
wave symmetry. However, the p-w&ve superconductivity 
can be stabilized in the present case owing to the orbital 
degeneracy. 

It should be noticed that the origin of trapezoidal peak 
of spin susceptibility around T point is clearly under- 
stood by this momentum dependence of wave function. 
Although the wave functions are not orthogonal between 
different pairs of hole pockets, the matrix elements be- 
tween them in calculating xil) are small. Therefore, in 
the zeroth order approximation, pairs of hole pockets are 
regarded to be decoupled from each other. Then, each 
hole pocket independently induces the trapezoidal peak 
of \(q) as in the two-dimensional electron gas model. 

Another point to stabilize the superconductivity is the 
disconnectivity of the e g -Fermi surface as discussed be- 
fore the discovery of Na x Co02 ■ yH 2 0. 4 °) Even in the 
anisotropic superconductivity such as p-wave or /i-wave 
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symmetry, the order parameter can take a same sign in 
each hole pocket, which stabilizes the superconductiv- 
ity induced by the ferromagnetic spin fluctuation. Note 
that the difficulty of the ferromagnetic spin-fluctuation- 
induced superconductivity (superfluidity) has been dis- 
cussed for 3 Hc. 54 ) This difficulty is removed by the topo- 
logical aspect of Fermi surface in case of Na x Co02-2/H 2 0. 

3.4 Momentum dependence of superconducting gap 

Next, we show the momentum dependence of order pa- 
rameter A(fe, mT) in Fig. 11. Although A e does not reach 
1 at T = 0.01 (Fig. 5), it is generally expected that the 
amplitude |A(fc,i7rT)| shows the momentum dependence 
of superconducting gap below T c and determines the low 
energy excitation. We note that even if the supercon- 
ducting instability is dominated by the e g -Fermi surface, 
the aig-Fermi surface also contributes to the low energy 
excitations observed by NMR 1/TiT, specific heat and 
magnetic field penetration depth. 

Fig. 11(a) shows the order parameter in the p-wave 
symmetry. We choose the Hund's rule coupling as Jh = 
U/3 where the p-wave superconductivity is stabilized. 
Among the two degenerate p x - and p y -states, only the 
Py-state is shown. Because of the discontinuity of the 
e g -Fermi surface, the order parameter is node-less on 
the eg-Fermi surface, while it has nodes on the ai g - 
Fermi surface. Since p x x ± p y y, p x y ± p y x or (p x ± ip y )z 
states are expected below T c , the superconducting gap 
becomes y/ A x (k) 2 + A y (fc) 2 , where A x (fc) and A y (k) are 
the order parameters for p x - and p y -states, respectively. 
In this case, the superconducting gap docs not vanish 
even on the ai g -Fermi surface. But, we find a remark- 
able anisotropy of the superconducting gap on the aig- 
Fermi surface which can explain the power-law behaviors 
of NMR 1/TxT and so on, like in the case of Sr 2 Ru0 4 . 55) 
However, we note that this is an accidental result. 

Fig. 11(b) shows the order parameter in the /i-wave 
symmetry. We choose the Hund's rule coupling as Jh = 
U/6 where the /i-wave superconductivity is most stable. 
We can see the clear six times alternation of the sign 
of order parameter. Also in this case, the e g -Fermi sur- 
face is node-less and ai g -Fermi surface has line nodes. As 
we showed before for the magnetic penetration depth, 15 - ) 
the combination of fully gaped e g -Fermi surface and line 
nodes on ai g -Fermi surface gives an intermediate temper- 
ature dependence between s-wave and anisotropic super- 
conductivity. 

In Fig. 11(c) we show the order parameter in the d xy - 
wave state which is stabilized when Jh is very small, 
Jh = U/12. The d xy ± iJ x 2_ y 2 state is expected below 
T c and both ai g -Fermi surface and e g -Fermi surface are 
node-less in this case. The exponential behaviors in many 
quantities are expected unless some accidental situation 
occurs as in the p-wave state. Our calculation does not 
support such an accidental situation in the d-wave sym- 
metry. 

It should be noticed that in all of the above cases we 
have shown, the amplitude of order parameter is large on 
the eg-Fermi surface, while it is small on the ai g -Fermi 
surface. This result is expected from the fact that the 
e g -Fermi surface is responsible for the pairing instability 
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Fig. 11. Momentum dependence of order parameter (a) in the 
Py-wave symmetry, (b) in the /i-wave symmetry and (c) in the 
dxy-wave symmetry. The parameters are chosen to be a = 0.8 
and n c = 0.31. 



as discussed in §3.3. This point will be illuminated more 
clearly in the next section. 

4. Reduced Models 

We have analyzed the possibility of unconventional su- 
perconductivity in Na x Co02 • J/H 2 on the basis of the 
three-orbital model. Because calculations for this model 
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need much computational time, a simplified model ap- 
propriate for studying the superconductivity is highly de- 
sired for a future development in the theoretical field. In 
this section, we try to find an appropriate model from the 
comparison to the three-orbital model. We show that the 
two-orbital model is satisfactory for this purpose, while 
the single-orbital model is not. The essential origin of the 
results in §3 will be clarified by these trials. 

4-1 Failure of single- orbital Hubbard model 

Thus far, we have stressed some essential roles of the 
orbital degeneracy. They arc illuminated by showing the 
failure of single-orbital model. Some authors have already 
studied single-orbital Hubbard models reproducing the 
LDA Fermi surface. 30 ' 31 -* In this paper, we try a single- 
orbital Hubbard model by keeping only the 7-band, 
i.e, the highest-energy eigenstates obtained in eq. (15). 
Hamiltonian is expressed in the following way. 



Hl = J2E 3 (k)cl s C ks + U^ 



(20) 



As has been shown in Fig. 1, the typical Fermi surface 
is reproduced in this model. Indeed, this is the minimal 
model describing the electron correlation in this material. 
However, as shown below, this model is inappropriate 
for the study of superconductivity because the results 
are qualitatively different from those in the multi-orbital 
model. 

We clarify the term "single-orbital Hubbard model" 
in order to avoid any confusion. In this paper, "single- 
orbital Hubbard model" suggests the single-band model 
including only the momentum independent interaction 
like eq. (20). As is shown later, we can construct a single- 
band model in which the roles of orbital degeneracy are 
appropriately represented in the momentum dependence 
of interaction term. Thus, we distinguish "single-orbital 
Hubbard model" from 'single-band model" . 
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Fig. 12. Phase diagram of the single-orbital Hubbard model. The 
qualitatively different results from Fig. 3 indicate the failure of 
this model. 



In Fig. 12, we show the phase diagram obtained by 
the SOP applied to the single-orbital Hubbard model 
(eq. (20)). In the wide region of parameter space, the 
d-wave and i-wave superconductivities are stabilized in- 
stead of p-wave and /i-wave states. The /i-wave super- 



conductivity competes with the d-wave one, but is stabi- 
lized only in a narrow region. The p-wave superconduc- 
tivity is not stabilized in the whole parameter range. 

This difference arises from the disregard of the momen- 
tum dependence of wave function which is represented by 
U(k). If we neglect the momentum dependence of U(k) 
in eq. (15), the three-orbital model is reduced to the 
single-orbital Hubbard model in eq. (20). The difference 
of stable pairing state is apparent if we check the spin 
susceptibility x{q)- I n the single-orbital Hubbard model, 
x(q) is similar to that obtained in Ref. 31 and we do not 
clearly see the ferromagnetic tendency (see also the dis- 
cussion in §3.3). As a result, the momentum dependence 
of the effective interaction is qualitatively different form 
that in the three-orbital model. 

This difference is partly improved by neglecting the 
aig-orbital like Ref. 30. Then, we obtain the nearly fer- 
romagnetic spin fluctuation and spin triplet supercon- 
ductivity. However, the coupling between different pairs 
of hole pocket Fermi surfaces (see Fig. 10(b)) is over- 
estimated, and therefore, the /x-wave state is stabilized 
much more than the p-wave state. This is not consis- 
tent with the results in §3. We wish to stress again that 
the characteristic nature of orbital in each hole pocket 
Fermi surface induces the nearly degeneracy between the 
p-wave and /i-wave states. This characteristic nature can 
not be taken into account in the single-orbital Hubbard 
model. 

4-2 Effective two-orbital model 

The results in the previous subsection show that the 
single-orbital Hubbard model is qualitatively inappro- 
priate for studying the superconductivity. The impor- 
tant factor to be taken into account is the orbital char- 
acter of quasi-particles on each Fermi surface. This is 
described by the momentum dependence of the unitary 
matrix U(k) in eq. (15). Considering these points, we 
propose a simplification of the three-orbital model in this 
subsection. The reduced model is an effective two-orbital 
model representing the e g -doublet. The simplification is 
performed by the following two steps. 

(i) The aig-orbital is simply ignored. 

(ii) The lower band below the Fermi level is ignored. 
The first step is justified because we find that the su- 
perconducting instability is dominated by the six hole 
pocket Fermi surfaces which mainly consist of the e g - 
orbitals. The second one is generally justified because 
the quasi-particles around the Fermi surface lead to the 
superconductivity. 

In order to perform the first step, we transform the 
basis of local orbitals. This is carried out by using the 
unitary transformation as, 

(4, >4 >4 ) = (4 ,ct ,4, )u h (21) 

v fe,l,s' fe,2,s' K,3,s J y fe,l,s' fc,2,s' K,3,s J " v ' 



V3 



4r 



1 

V2 




(22) 



The interaction term Hi in the Hamiltonian H3 is invari- 
ant for this unitary transformation owing to the relations 
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U = U' + 2 Jh and Jh = J. The non-interacting term is 
transformed as, 



k,s 

H'{k) = u}H{k)U h 



(23) 



(24) 



The first step is performed by dropping the creation (an- 
nihilation) operator d^ i (d k 1 g ) which corresponds to 

the aig-orbital. As a result, the three-orbital model is 
reduced to the following two-orbital model. 

2 

k,s 1 a=1 



+U' £ £ rii^n^b - J H £ £(2S ii0 S i)(l 



(25) 



a=tb 



H'(k)i+ij+i and two component vector aj^ 
(rfj, o , dj^ ). Then, the Green function is described by 



Here, we have introduced a 2 x 2 matrix h(k)i j — 

a 2x2 matrix as G(k) = (iw n l — /i(fc)) _1 , whose elements 
are expressed as 

2 

Gij(k) = J2v ia (k)v ja (k)G a (k). (26) 

a=l 

Here, Vi a (k) are components of the unitary matrix V^{k) 
which diagonalizes the matrix h{k) 



V\k)h{k)V{k) 



ei(fc) 
e 2 (k) 



(27) 



with ei(fe) < e 2 (fc). The diagonalized Green function is 
obtained as G a (k) = {uJn ^ e ^ ky 

We show the dispersion relation ei(fc) and e 2 (fc) in 
Fig. 13. Apparently the band structure around the e g - 
Fermi surface is unchanged by this simplification, while 
the aig-Fermi surface vanishes. 
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Fig. 13. Dispersion relation in the two-orbital model (solid lines). 
The parameters arc chosen to be a = 0.8 and n c = 0.36. We have 
shown the dispersion relation of the ai g -orbital which is obtained 
as fT (fe)n — fi (dashed line). 



The second step is performed by ignoring the lower en- 
ergy band, ei(fc). Then, the Green function is obtained 
as, Gij(k) — Vi2(k)vj2{k)G2(k). Owing to this proce- 
dure, the calculation becomes equivalent to that for a 
single band Hamiltonian with momentum-dependent in- 
teraction, 



E 

k,s 



e 2 (fc)ct 



k ' k 



+ £ S{q,k' k / k . f k > A c kA 
q,k',k 



+ S'(q,k',k)cl 



k. 



-fe',<r k',a c k 



(28) 



q,k',k,a 



The momentum dependent factors S(q, k' , k) and 
S'(q, k! , k) are expressed by the Coulomb interactions U, 
U' , Jh and J and the wave function Vi2(k). If we neglect 
the momentum dependence of unitary matrix V\k), the 
factor S(q,k',k) becomes U and S'(q,k',k) — 0. Then, 
the model is exactly reduced to the single-orbital Hub- 
bard model described by eq. (20) with use of e 2 (k) in- 
stead of E 3 (k). We have discussed in §4.1 that this single- 
orbital Hubbard model is not appropriate. On the other 
hand, the Hamiltonian Hl is appropriate because the 
roles of orbital degeneracy are taken into account in the 
momentum dependence of interaction. 
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Fig. 14. Eigenvalues of Eliashbcrg equation obtained in the effec- 
tive two-orbital model. We choose the parameters a = 0.8 and 
J H = (7/4. 



We find that the results for the superconductivity are 
almost the same between the Hamiltonian H 2 and Hl- 
In Fig. 14 we show the n c -dependence of eigenvalues of 
Eliashberg equation for the simplified model, Hl- We see 
that the increase of eigenvalues with n c is steeper than 
that in Fig. 4. This is mainly owing to the increase of 
DOS. However, the relation between each pairing sym- 
metry closely resembles. For example, the p-wave super- 
conductivity is stable around n c — 0.2, while the /i-wave 
superconductivity is realized for larger values of n c . The 
eigenvalue for the spin singlet <i-wave superconductivity 
is far below that for the spin triplet one. These results 
mean that the effective two-orbital model described by 
cq. (25) or eq. (28) appropriately reproduces the results 
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in the three-orbital model. The fact that the step (1) 
is appropriate clearly means that the superconductivity 
is basically led by the e g -Fermi surface. The ai g -Fcrmi 
surface plays only a secondary role. 

Note that the eigenvalue of Eliashberg equation de- 
creases owing to the step (1), mainly because the DOS in 
the e g -Fermi surface decreases. We have confirmed that 
the step (2) slightly enhances the spin triplet supercon- 
ductivity. 

5. Effects of Vertex Corrections in a Two- 
Orbital Model 

In this section, we study the effects of vertex correc- 
tions. Although it is desirable to study these effects in 
the three-orbital model, we use the effective two-orbital 
model whose validity has been demonstrated in §4.2, be- 
cause of numerical difficulties. Generally speaking, the 
higher order terms may play an important role for the 
superconducting instability, since it is considered that 
most of unconventional superconductors are in the inter- 
mediate coupling region. For example, vertex correction 
which is not included in the RPA plays an important 
role to stabilize the spin triplet pairing in Sr 2 Ru04. 36 ) 
Therefore, it is an important issue to investigate the role 
of higher order corrections in the present model. 
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Fig. 15. Diagrammatic representation of the third order terms 
in the effective interaction, (a-f) correspond to the spin singlet 
channel or spin triplet channel with d-vector d \\ z. (c'-f) corre- 
spond to the spin triplet channel with d-vector d _l_ z. 



We apply the third order perturbation theory (TOP) 
and its renormalized version to the Hamiltonian H2 
(cq. (25)). We adopt this model instead of more sim- 
plified model Hi (eq. (28)) because the computational 
time is hardly reduced by the second step (ii) in §4.2. 
The parameter is chosen to be Jh = U/3, where the in- 
teraction between electrons with same spin vanishes and 
thus the number of diagrams is much reduced. As dis- 
cussed in §3.2, this region will be relevant rather than 
the region where the Hund's rule coupling is small. 

Fig. 15 shows the diagrammatic representation of third 
order terms in the effective interaction. Figs. 15(a) and 



(b) are classified into the RPA terms and others are the 
vertex corrections. The present theory is invariant for 
the rotation of spin, since we do not take account of the 
spin-orbit interaction. Therefore, the result on the spin 
triplet pairing does not depend on the direction of d- 
vector. Note that two RPA terms cancel each other in 
case of the spin triplet pairing with d\\ z. 



p-wave (SOP) 

• • p-wave (TOP) 

^--^fj-waveCTOP) 
A---Af 2 -wave(TOP) 
* i-wave (TOP) 




Fig. 16. Eigenvalues of Fliashbcrg equation in the third order 
perturbation theory. The thick solid line shows the maximum 
eigenvalue in the second order perturbation theory, which is clas- 
sified into the p-wave symmetry. We do not show the eigenvalue 
in the c(-wave symmetry because the tendency to superconduc- 
tivity is very weak. We fix the parameters a = 0.6, n c = 0.35 
and J H = U/3. 



We numerically solve the Eliashberg equation within 
the TOP and show the eigenvalues in Fig. 16. We see 
that the p- and /2-wave superconductivity are signifi- 
cantly stabilized for U > 4, while the /i-wave and spin 
singlet pairings are unfavored. However, as discussed be- 
low, we find that these results in the intermediate cou- 
pling region are fictitious. Within the third order terms in 
Fig. 15, dominant contributions for triplet channel come 
from the terms represented in Figs. 15(e') and (f), which 
include a particle-particle ladder. In contrast, the terms 
represented in Figs. 15(c') and (d') with a particle-hole 
ladder are negligible. As is well known in the Kanamori 
theory on the metallic ferromagnetism, 56 ) the particle- 
particle ladder diagrams generally induce the screening 
of interaction as U — ► U(q) = U/ (1 + U<p(q)) where <j)(q) 
is obtained by the particle-particle ladder diagram. If 
g-dependence of U(q) is not important, this scattering 
process is incorporated by the renormalized coupling con- 
stant U. In the above TOP calculation, only the lowest 
order term in the Kanamori-type correction was taken 
into account. Therefore, it is reasonable to think that the 
contributions from Figs. 15(e') and (f ') can be suppressed 
if we include the higher order perturbation terms. 

In order to investigate this possibility, we perform a 
calculation of a renormalized TOP (RTOP), as shown in 
Fig. 17. The particle-particle ladder in Figs. 15(e') and 
(f ') are replaced by the T-matrix shown in Fig. 17(a). As 
a result, the infinite order terms representing the screen- 
ing effect are taken into account as in the Kanamori 
theory. By using the diagrams in Figs. 2, 15(c',d') and 
17(e",f"), we estimate the effective interaction and solve 
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Fig. 17. (a) Renormalization of the particle-particle ladder dia- 
gram. In (e") and (f"), the renormalized particle-particle ladder 
is used instead of bare ladder. In the RTOP, we take account of 
the terms (e") and (f") instead of (e') and (f) in Fig. 15. 



the Fliashberg equation. The obtained eigenvalues are 
shown in Fig. 18. It is apparent that the results of naive 
TOP is significantly altered by the renormalization and 
that the correction to the SOP is small. In particular, 
the p-wave superconductivity is slightly stable over the 
/i-wave superconductivity. The nearly degeneracy be- 
tween these states is also reproduced. The order param- 
eter in each pairing symmetry is very similar to Fig. 11, 
although that in the naive TOP is remarkably different. 
We see that the eigenvalues are slightly reduced from the 
SOP, however the {/-dependence is almost unchanged. 
These results are naturally interpreted if we consider that 
the vertex corrections basically work as a screening ef- 
fect. Then, the second order perturbation theory is justi- 
fied by regarding the interactions to be the renormalized 
ones. 



1 
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'0.4 
0.2 



p-wave (RTOP) 
<y-0t , -wave (RTOP) 
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Fig. 18. Eigenvalues of Eliashberg equation in the renormalized 
third order perturbation theory for p-wave (circles) and /i-wave 
(diamonds) symmetry. Note that the eigenvalue in the /2-wave 
symmetry is very small. The thick solid and dashed lines show 
the eigenvalues in the second order perturbation theory for the 
p-wave and /i-wave symmetry, respectively. The parameters are 
the same as in Fig. 16. 



Let us compare the present results to the case of high- 
T c cuprates and Sr 2 Ru04. For high-T c cuprates, the d- 
wave superconductivity is basically induced by the RPA 
terms and the vertex correction due to the particle- 
particle ladder diagrams effectively reduces the coupling 
constant. 35, 57 -* Therefore, the situation is very similar to 
the present case, although there is a difference of sin- 
glet and triplet pairing. On the other hand, in case of 
Sr 2 Ru04, the effective interaction derived from the RPA 
terms has very weak momentum dependence, which does 
not work for the anisotropic pairing. However, the in- 
dependence of particle-particle ladder in TOP favors the 
spin triplet superconductivity. 36 -* Then, the naive discus- 
sion on the screening effect can not be applied. It has 
been confirmed that the qualitative results of TOP ap- 
plied to Sr2Ru04 are not altered even when the renor- 
malization of particle-particle ladder is taken into ac- 
count. 58 ) Thus, the basic mechanism of possible spin 
triplet superconductivity in Na x Co02 • J/H 2 is quali- 
tatively different from that in Sr 2 Ru04. 

6. Discussions 

In this paper, we have investigated the multi-orbital 
model for Na x Co0 2 • yH 2 on the basis of the perturba- 
tion theory. The obtained results indicate a possibility of 
spin triplet superconductivity in this material, although 
the G?-wave superconductivity is also stabilized in a part 
of parameter space. There are two candidates of spin 
triplet pairing; p-wave and /-wave superconductivity are 
nearly degenerate. 

Although the spin triplet superconductivity is one 
of the most interesting issues in the condensed matter 
physics, the microscopic theory remains in the develop- 
ing stage. This is mainly owing to very few (i-electron 
materials showing the spin triplet superconductivity. Al- 
though we see many candidates in the heavy fermion ma- 
terials, the theoretical treatment is generally difficult for 
/-electron systems. Therefore, a discovery of spin triplet 
superconductor in transition metal oxides will lead to an 
important development in the microscopic understand- 
ings. 

Probably, most established spin triplet superconduc- 
tor in (i-electron systems is Sr 2 Ru04. 34 - ) Therefore, we 
have provided detailed discussions on the comparison be- 
tween Sr 2 Ru04 and Na x Co0 2 • yH 2 0. According to the 
results in this paper, Na x Co0 2 • yH 2 provides a quali- 
tatively different example from Sr 2 Ru04 in the following 
two points. 

First, the RPA terms give rise to the dominant scatter- 
ing process leading to the spin triplet pairing. The spin 
excitation is clearly ferromagnetic and favorable for the 
spin triplet pairing. This is in sharp contrast to the case 
of Sr 2 Ru04 where the vertex corrections are essential for 
the p-wave pairing. In case of Na x Co0 2 • J/H 2 0, the ver- 
tex corrections induce only the screening effect which is 
not important for the qualitative results. While the ferro- 
magnetic spin-fluctuation-induced spin triplet supercon- 
ductivity has been discussed from early years, the cor- 
responding superconductivity has not been established 
until now. We expect that Na x Co0 2 • yil 2 will be a 
first example realizing this mechanism. 
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Second, the orbital degeneracy plays an essential role 
in case of Na x Co02 • 2/H2O. The conduction band in 
Na x Co02 • yH 2 as well as that in Sr 2 Ru04 arc ba- 
sically described by three t 2g -orbitals. Although the 
single-orbital Hubbard model is an appropriate model 
for describing the pairing mechanism of S^RuG^, 36 ** 
such a simplification is qualitatively inappropriate for 
Na x Co02 • yH 2 0. The success of single-orbital Hubbard 
model for Sr 2 Ru04 is due to the electronic structure 
where the 7-band is basically described by the local d Ky - 
orbital. The failure for Na x Co02-2/H 2 is due to the fact 
that the e g -Fermi surface can not be described by any in- 
dividual local orbital. In other words, the hybridization 
term in the unperturbed Hamiltonian is large in case 
of Na x Co02 • J/H2O, while it is negligible in Sr2Ru04 
owing to the particular crystal symmetry. In this sense, 
Na x Co02 • J/H 2 will be a more typical example of 
the multi-orbital superconductor. Then, the momentum 
dependence of the wave function of quasi-particles es- 
sentially affects the effective interaction leading to the 
Cooper pairing. 

We have pointed out that the reduced two-orbital 
model is appropriate, instead of the failure of single- 
orbital model. This is because the Fermi surface in 
Na x CoC-2 • yH 2 can be classified according to the local 
orbitals. Then, the superconductivity is basically trig- 
gered by the e g -Fermi surface. Since a portion of ai g - 
orbital in the e g -Fermi surface is less than 5%, this orbital 
is safely ignored. This situation is similar to the case of 
Sr 2 RuC>4. However, the orbital degeneracy in e g -doublet 
can not be ignored in case of Na x Co02 • yH 2 0. 

From the above comparisons, we obtain the following 
empirical rules. 

(1) When the RPA-terms are favorable for the 
anisotropic superconductivity, the non-RPA terms are 
not qualitatively important, and vice versa. 

(2) When a part of Fermi surface is described by a few 
local orbitals, the simplification of microscopic model is 
possible. 

In particular, the second rule will be helpful for a fu- 
ture development of microscopic understanding on the 
multi-band superconductors. For example, several Fermi 
surfaces appear in heavy fermion materials. This fact as 
well as the 14-fold degeneracy in /-shell make the micro- 
scopic treatment difficult. However, it will be possible to 
obtain a simplified model by identifying the microscopic 
character of each Fermi surface. 

Thus far, we have discussed the superconductivity in 
Na x Co02 • J/H2O induced by the electron-electron cor- 
relation and highlighted the possibility of spin triplet 
pairing. However, any clear experimental evidence for 
the symmetry of superconductivity has not been ob- 
tained up to now. Instead, we see some experimental 
observations which restrict the pairing state. For exam- 
ple, the absence of (or very small) coherence peak in 
NMR l/TiT, 10 ~ 13 ) power-law temperature dependence 
of l/TiT 12 < 13 ) and specific heat, 17 " 19 ) NMR Knight shift 
below T c i ' 11 . 59 ' 60 ) an( j time-reversal symmetry ob- 
served in /iSR 14 ' should be cited, although a part of 
them are controversial. As for the results in this paper, 
spin triplet p- or /i-wave superconductivity is consistent 
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with the absence of coherence peak and with the power- 
law behaviors below T c . In both cases, the (quasi-)line 
nodes appear in the ai g -Fermi surface. In case of the 
p-wave pairing, the time-reversal-symmetry observed in 
[iSR indicates a <i-vector parallel to the plane, namely 
d = p x x ± p y y or d = p K y ± p y x. This direction of 
d- vector is consistent with the recent measurements of 
NMR Knight shift under the parallel field 11 ' 59 ' 60 ) as well 
as macroscopic H C 2, 7 ' 61 ) if we assume that the d- vector is 
strongly fixed against the applied magnetic field. We note 
that the qualitatively different result has been obtained 
in the NMR Knight shift, 10 ) which is consistent with this 
pairing state if the <i-vector is weakly fixed against the 
magnetic field. 

Although we have shown that the d- vector in Sr 2 Ru04 
is very weakly fixed against the magnetic field, 37 ) this 
is partly owing to the particular electronic structure of 
Sr 2 RuC"4. Therefore, we expect that the anisotropy of 
d-vector is larger for Na x Co02 • J/H2O. The symmetry 
breaking interaction leading to the anisotropy arises from 
the second order term with respect to the spin-orbit in- 
teraction for Sr2RuC>4, while it arises from the first order 
term in case of Na x Co02 • yH 2 0. Therefore, it is possible 
that the cZ-vector is strongly fixed against the magnetic 
field in case of Na x Co02 • 2/H2O. Quantitative estima- 
tions for the anisotropy will be one of the interesting 
future issues. 

On the other hand, the possibility of spin singlet su- 
perconductivity has not been denied up to now. Then, 
the absence of time-reversal symmetry breaking will be 
a issue to be resolved for d-wave pairing because the 
d x 2_ y 2 ± id xy state is expected so as to gain the conden- 
sation energy. The local distortion of triangular lattice 
or the feedback effect will be a candidate of the reso- 
lution. It seems that the i-wave superconductivity 62 ) is 
consistent with the present experimental results except 
for the very weak impurity effects. 63 ) However, the mi- 
croscopic mechanism leading to the pairing with T c = 5K 
will be difficult for such a high angular momentum state. 
In our study, we have not found the stable i-wave state. 
Although the observed impurity effect seems to support 
the s-wave pairing which is robust for the disorder, very 
short quasi-particle life time or significant anisotropy in 
the gap function has to be assumed for the absence of co- 
herence peak in \jT\T . We consider that further vigorous 
investigations are highly desired for the identification of 
pairing state in Na x Co02 • J/H 2 0. 
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